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> ! Abstract 
m 

^Sj ' We study a double inflation model (a hybrid inflation + a new inflation) 

' in supergravity and discuss the formation of primordial black holes (PBHs) 

. with mass ~ 10"'^'^ — 10^ M0. We find that in a wide range of parameter 
space, we obtain PBHs which amount to ~ 1, i.e., PBH dark matter. Also, 

Qh! we find a set of inflation parameters which produces PBHs evaporating now. 

pi ^ Those PBHs may be responsible for antiproton fluxes observed by the BESS 

^ ' experiment. 
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J3 



I. INTRODUCTION 

In the framework of supergravity the reheating temperature of inflation should be low 
enough to avoid overproduction of gravitinos |]T],0]. The new inflation model generally 
predicts a very low reheating temperature and hence it is the most attractive among many 
inflation models |^ . However, the new inflation suffers from a fine-tuning problem about the 
initial condition; i.e., for a successful new inflation, the initial value of the inflaton should 
be very close to the local maximum of the potential in a large region whose size is much 
longer than the horizon of the universe. 

A framework of a double inflation ||^ was proposed to solve the initial value problem of 
the new inflation modelQ. It was shown that the above serious problem is naturally solved by 
supergravity effects if there exists a preinflation (e.g., hybrid inflation 0) with a sufficiently 
large Hubble parameter before the new inflation p[. 



^Different models of double inflation were studied by various authors 
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In this double inflation model, density fluctuations produced by both inflations are cos- 
mologically relevant if the e-fold number of the new inflation is smaller than ~ 60 (the 
total e-fold number ~ 60 is required to solve flatness and horizon problems in the standard 
big bang cosmology In this case, the preinflation should account for the density fluc- 
tuations on large cosmological scales [including the Cosmic Background Explorer (COBE) 
scales] while the new inflation produces density fluctuations on smaller scales. Although the 
amplitude of the fluctuations on large scales should be normalized to the COBE data [Q, 
fluctuations on small scales are free from the COBE normalization and can have arbitrary 
power matched to observations. In Ref. [rU|, a cosmological implication of the double infla- 



tion for the large-scale structure formation was discussed. In this paper, we study primordial 
black hole (PBH) formation in the double inflation model. In Refs. []TTl , p!2| , the production 
of black hole MACHO was investigated in the double inflation model for a special caseQ. 
Here, we consider a wide range of parameter space where PBHs are formed^. In particular, 
we show that the double inflation creates small PBHs evaporating now if those PBHs are 
produced during matter-dominated (MD) era, i.e., before the end of reheating process after 
the new inflation. We stress that these evaporating PBHs may account for antiproton fluxes 



observed by the BESS experiment [14 



Throughout this paper the gravitational scale (~ 2.4 x 10 GeV) is taken to be unity 



II. BLACK HOLE FORMATION 



In a radiation-dominated (RD) universe, PBHs are formed if the density fluctuations 5 
at horizon crossing satisfy a condition 1/3 < 6 < 1 \TE,W\, where 6 is the over density at 



the horizon scale. Masses of the black holes M are roughly equal to the horizon mass, 

R / T / a \ ^1/2 

M ^ 47r. - ^ 0.066Me — - |^ , (1) 



p ■ VGeV/ V50 

where p, T, and g^, are the total cosmic density, temperature, and statistical degrees of 
freedom at the horizon crossing, respectively. 

The horizon length at the black hole formation epoch [T = T^) corresponds to the scale 

in the present universe given by 

...|^,-,.,.0.0e4pc(^)-'(|)-, 



^Different models for the PBH formation have been studied in Ref. [13|. 

^In this paper we investigate the PBHs with mass ~ 10"^'' — IO^Mq. The upper bound of 1O^M0 
is not rigorous but much heavier black holes can be produced in the present model if we take 
the appropriate model parameters. However, the mass of the black holes should be less than the 
galactic mass ~ lO^^M©, otherwise the power spectrum conflicts the observations ( e.g. distribution 
of galaxies ). The lower bound ~ 10~^'^Mq comes from the requirement that e-fold number of the 
new inflation should be larger than 0. 
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where Tq is the temperature of the present universe. The comoving wave number corre- 
sponding to this length scale, = 27i/L^, is 



...1.0xlO=Mpc-.(|)'"(|^). 
Thus, we can write the PBH mass as a function of comoving wave number as 



(3) 



M..6.4xlO"Me(|)"75i|^)". (4) 



The mass fraction /5*(= Pbh/p) of PBHs of mass is given by [jT6 



/^*(^*) = / , ^=T77Texp (-TT-^TTT ) - ^(M.) exp (-—^-—] , (5) 



1/3 v^(t(M,) " V 2a2(M,) y - V - \^ 18a2(M0 

where cr(M*) is the mass variance at the horizon crossing. Notice that the mass fraction 
/3*(M^=) drops off sharply as cr(M) decreases. The density of the black holes of mass M^=, 
Pbh{M^), is given by 

'-^^ ^ ^-MM.)T., (6) 
s 4 

where s is the entropy density. Since Pbh/s is constant at T < T*, we estimate the density 
parameter Qbh{M^) of the black holes in the present universe as 

nBH{M,)h' ~ 2.1 X 10V*(M,) (^^) , (7) 

where h is the present Hubble constant in units of 100 km/sec/Mpc. We write it as a 
function of PBH mass or PBH scale as 

QBH{M^)h' ^ 5.4 X 10^/3.(M,) (|)"'^' (^) ^'^ , (8) 



or 



^Bn{M:)h^ - 2.1/?.(M.) (I)"'' . (9) 

As for the mass of the PBHs produced during the RD era, we have a lower limit Mr. 
The mass of PBH produced at the reheating epoch is given by [see Eq.(|l|)] 

M^^^mmJ^\ (10) 



GeV/ V50 

where Tr is the reheating temperature. As seen later is less than lO^GeV in our inflation 
model, and hence Mr is larger than ~ 10^^'^ M©. Thus, the PBHs lighter than M/j should 
be produced during the MD era. 

In a MD universe, a relation between the comoving scale and the horizon mass M^, is 
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L*=LrU^) , (11) 



R 



where Lji is the comoving scale of the horizon at the reheating epoch. Thus, the mass 
of the PBH produced during the MD epoch is given by 



^ 6.3 X 10^^ Mf. 







VGeV/ iMpc 



We see that small PBHs of mass, ~ 10 ^^Mq for example, would be produced for 
Tr ~ lO^GeV and /c^, ~ lO^^Mpc"^. The condition for the PBH formation in a MD universe 



is discussed in Ref. |17|, where the mass fraction of PBHs of mass is estimated as [17 



(3,{M,) ~ 2 X 10-V(M,)i='/2_ (^3) 

Notice that the mass fraction has a weeker dependence on a than in the RD case [see Eq.(^)] 
During the MD era, the mass fraction (3^ stays constant and hence the density of the 
black holes of mass M=k, Pbh{M^), is given by 

^ ^-^{M^)Tn. (14) 
s 4 

We can write the present density parameter Qbh{M^,) for the black holes of mass as 

QBHiM,)h^ ~ 2.1 X 10V*(M,) (^^) . (15) 



III. DOUBLE INFLATION MODEL 

We adopt a double inflation model proposed in Ref. The model consists of two 



inflationary stages; the first one is called preinfiation and we take a hybrid inflation [|T8 
(see also Ref. [1^) as the preinfiation. We also assume that the second inflationary stage is 
realized by a new inflation model [^ and its e-fold number is smaller than ~ 60. Thus, the 



density fluctuations on large scales are produced during the preinfiation and their amplitude 
should be normalized to the COBE data [Q. On the other hand, the new inflation produces 
fluctuations on small scales. Since the amplitude of small scale fluctuations is free from the 
COBE normalization, we expect that the new inflation can produce large density fluctuations 
enough to form PBHs. We choose the predicted power spectrum to be almost scale invariant 
{us ^ 1) on large cosmological scales which is favored for the structure formation of the 
universe [^. On the other hand, the new inflation gives the power spectrum which has 
large amplitude and shallow slope (n^ < 1) on small scales. Thus, this power spectrum has 
a large and sharp peak on the scale corresponding to a turning epoch from the preinfiation 
to the new inflation, and we expect that PBHs are produced at that scale. 



As for the detailed argument of the dynamics of our model, see Refs. [10-12 
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A. Preinflation 



First, let us briefly discuss a hybrid inflation model ||18|. The hybrid inflation model 
contains two kinds of superfields: one is S{x,9) and the others are a pair of \E'(x, 6^) and 
^{x,6). Here 6 is the Grassmann number denoting superspace. The model is based on the 
U(1)r symmetry under which S{e) e2*"^(^e-*") and ^(^)ff(e) ^(ee-*")^(^e-^°). The 
superpotential is given by [|T8| 



W{S,^,<il) = -fi'^S + XS^^. (16) 

The i?-invariant Kahler potential is given by 

K{S, ^, ^) = \S\^ + + l^rp + . . . , (17) 

where the ellipsis denotes higher-order terms which we neglect in the present analysis for 
simplicity. We gauge the U(l) phase rotation:\E' e^^"^ and 'I' e"*''^'. To satisfy the 
D-term flatness condition we take always \E' = ^ in our analysis. 

We define A'^cobe as the e-fold number corresponding to the COBE scale and the COBE 
normalization leads to a condition for the inflaton potential, 



Y3/2 



v 



^ 5.3 X 10"^ (18) 

A^COBE 



where V is the inflaton potential obtained from Eqs.(|T^) and (p^Tf ). In the hybrid inflation 
model, density fluctuations are almost scale invariant; 



-^COBE 

jVcobe 



^ 1, (19) 



where spectral index for a power spectrum of density fluctuations. 

B. New inflation 

Now, we consider a new inflation model. We adopt an inflation model proposed in 
Ref. The inflaton superfield (t){x,6) is assumed to have an R charge 2/{n + 1) and 

U(l)/i; is dynamically broken down to a discrete Z2nR at a scale f , which generates an 
effective superpotential [pO|,p|] , 



W{<P)=v^-^r^\ (20) 
The i?-invariant effective Kahler potential is given by 

m,x) = w + 1w+---. (21) 

where k is a constant of order 1. We require that supersymmetry breaking effects make the 
potential energy at a vacuum vanish, and we have a relation between v and the gravitino 



mass m3/2 as (for details, see Ref. P0| ) 



5 



The inflaton (f){x) (the scalar component of (j){x, 6)) has a mass m,^ in the vacuum with 

~ (23) 

The inflaton cj) may decay into ordinary particles through gravitationally suppressed inter- 
actions, which yields reheating temperature Tr given byQ 

Tr ~ O.lmJ/^ ~ 2.4 X 10^^GeVr2=^/2|^|2/2n|^|3{i~i/n) ^34) 
< lO^GeV for mg/a < ITeV, n > 3. (25) 

An important point on the above density fluctuations is that it results in a tilted spectrum 
with spectral index n^ew given by (see Refs. p|,pO|) 



(26) 



C. Initial value and fluctuations of the inflaton Lp 

The crucial point observed in Ref. is that the preinflation sets dynamically the initial 
condition for the new inflation. We identify the inflaton field Lp{x)/\/2 with the real part 
of the field 0(x). It gets an effective mass meg ~ /x^ during the preinflation 0. Thus, this 
inflaton ip tends to the potential minimum. 



V2 



V 



V - . (27) 



vA 

Notice that (pmin deviates from zero due to the presence of a linear term v'^^'^Sip (see Ref. 0). 
Thus, at the end of the preinflation the ip settles down to this v^min- 

After the preinflation, the universe becomes MD because of the oscillation of the inflaton 
for preinflation. During the MD era between the two inflations, the energy density scales 
as oc a~^, and the new inflaton oscillates around = with its amplitude decreasing as 
oc a"^/^. Since the scale factor increases by a factor (fx/v)^^'^ during this era, the mean initial 
value (pb of (p at the beginning of the new inflation is written as 

Therefore, the amplitude of fluctuations with comoving wavelength corresponding to the 
horizon scale at the beginning of the new inflation is given by 



The decay rate of the inflaton (p is discussed in Ref. 
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In V 



2/3' 



V 



-3/2 



-) 

V J 



4/3' 



-3/4 



3i/427r U 



(29) 



where ifpre is the Hubble parameter during the hybrid inflation, iJpj.^ ~ /i^/3. The fluctu- 
ations given by Eq. (|29| ) are a little less than newly induced fluctuations at the beginning 
of the new inflation [difnew — v'^ /{2n^/3)]. Moreover, the fluctuations produced during the 
preinflation are more suppressed for smaller wavelength. Thus, we assume that the fluctua- 
tions of induced in the preinflation are negligible compared with fluctuations produced by 
the new inflation. As mentioned before, the new inflation gives the tilted spectrum on small 
scales [see Eq. (p6D] and hence the fluctuations at the scale corresponding to the beginning 
of the new inflation is dominant. 

Now let us estimate e-fold number which corresponds to our current horizon. The e-fold 
number is given by |2^ 



k end r reh 

where \4 is a potential energy when a given scale k leaves the horizon. Vend that when the 
inflation ends, and preh energy density at the time of reheating. 

We take Vk — Vend, and p^^h — ^ ^ ^reh- For k = aoHo (i.e., the present horizon 
scale), we have 

A^tot^67.1+ (l--)\nv + l\nn + ^\ng. (31) 
V3 nJ 2 2n 

In estimating A^^cobe we must take into account the fact that the fluctuations induced at e- 
fold number less than (2/3) \n{fi/v) before the end of the hybrid inflation reenter the horizon 
before the new inflation starts. Such fluctuations are cosmologically irrelevant since the new 



inflation produce much larger fluctuations |TD[. Thus, A'cobe is given by 

2 /i 

iVcOBE = Atot - A'new + 77 ^ - 

3 V 

-67.1+ (^^-^ynv + ^\nn+^lng-Nne^+'^\n^. (32) 
The COBE normalization in Eq. ([T8| ) should be imposed by using this A'cobe- 



D. Numerical Results 

We estimate density fluctuations in the double inflation by calculating evolution of ip 
and a numerically. For given parameters k and A, we obtain the break scale kb and the 
amplitude of density fluctuations produced at the beginning of new inflation 5b- Here, k^^ is 
the comoving scale corresponding to the Hubble radius at the beginning of the new inflation 
(a turning epoch). We can understand the qualitative dependence of {kb,Sb) on (k. A) as 
follows: When k is large, the slope of the potential for the new inflation is too steep, and 
the new inflation cannot last for a long time. Therefore, the break occurs at smaller scales. 
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As for Sb, we can see from Eq.([T8D that as A gets larger, fi also gets large. In addition, we 
can show that 

oc ~ , (33) 

pi, K n 

for a given v (see Ref. |T^). Thus, we have larger 6b for larger A. 

IV. MASS VARIANCE AND DENSITY FLUCTUATIONS 

Our double inflation model predicts the amplitude of density fluctuations 6b as a function 
of inflation parameters. On the other hand, the black hole abundance is expressed as a 
function of mass variance a at the time when the corresponding scale enters the horizon. 
Therefore, when we compare the observations with the prediction of our model, we need a 
relation between the mass variance a and the fluctuations 6b- 

For the power spectrum with the break scale k^^ which enters the horizon during the 
RD epoch, we have a relation between the mass variance and the amplitude of fluctuations 

6b ~ (Tb/0.65, (34) 

(the numerical factor depends on the tilted spectral index Ug, and within the parameter 
range we consider, this factor lies between 0.62 ~ 0.67.). For the power spectrum with the 
break scale fc^^ which enters the horizon during the MD epoch, we have a relation between 
the mass variance and the amplitude of fluctuations as 

6b ^ crb/2.3, (35) 

(again, the numerical factor depends on the tilted spectral index Ug, and within the param- 
eter range we consider, this factor lies between 2.1 ~ 2.8.). 

First, let us consider PBH dark matter with Qbh ~ 1 which are produced during the RD 
epoch. Since the density fluctuations at the break scale is dominant, and the mass fraction 

has a sharp peak at that scale, only the PBHs of mass corresponding to the break scale 
are formed. For PBHs produced during the RD epoch (after reheating process), we have, 
from Eq. (|), 

n.„.-5.4xl0'ft(|)-"'g)"^ (36) 

For example, if we require that the black holes with mass ~ Mq (=MACHOs) be dark 
matter in the present universe, i.e. VLsnh^ ~ 0.25, we obtain ~ 5 x 10~^, and from Eq.(^ 
we obtain 



^In Refs. [^:12], an incorrect relation 5b — t^/G was used. 
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ct(M0) ~ 0.06. (37) 

From Eq.(^), we see that the break amphtude is 6b — 0.06/0.65 ~ 0.092, and we find from 
Eq.(H) that the break scale is 2.5 x lO^Mpc"^. In Fig.|l|, we plot the numerical results of our 
double infiation model for n = 4, g = 1, and v = 10^^. We see a wide range of parameter 
space which may account for DM {Qbh — 0.1 — 1). 

Next, let us consider another interesting mass range of PBHs which are evaporating now 
(Mevap ~ 3 X 1O~^^M0). The PBHs of such light mass are produced during the MD epoch 
and from Eq.(|T5|) we find 

QBH{M,)h' ^ 4.2 X 10V(M,)i3/2 (^Tn_^ ^^^^ 



It has been reported, recently, that the BESS experiment |Tj] has observed antiproton fiuxes, 
which may be explained by the evaporation of PBHs if VlBuh^ ~ 2 x 10~^ []T^. In order to 
explain the BESS result by evaporating PBHs we need 

/ \ —2/13 

a(Mevap) ^ 4.4 X 10-=^ j . (39) 

From Eq.(l^), we estimate the required fraction of the evaporating PBHs as 

C^(Me,ap) ^ 9.3 X 10-6^-3/13 1^ I -3/13n^-6(l-l/n)/13_ (40) 

Since the mass variance cr(M) scales as a{M) oc M'^-'^-"'^/^ |]ig| during the MD epoch, we 



obtain the mass variance at the break scale as 

(l-ns)/6 



9.3 X 10-6^-3/13|^|-3/13n^-6(l-l/n)/13 ^ (4^) 

\ -'"evap / 



and the amplitude is 5^ ~ cr;,/2.3 and ~ 1 — 2k. In Fig.^ we plot an example of the 
numerical results of our double infiation model^ for n = 3, g = 10"*^, and v = 10~^. As 
shown in the figure, we have a set of infiation parameters {k, A) which may account for the 
BESS experiment. 

V. CONCLUSIONS AND DISCUSSIONS 

In this paper we have studied the formation of PBHs by taking a double infiation model 
in supergravity. We have shown that in a wide range of parameter space PBHs are produced 
of various masses. These PBHs are interesting since, for example, they may be identified 
with MACHOs (M ~ Mq) in the halo of our galaxy. Or, they may be PBHs which are just 
evaporating now (M ~ IO'^^Mq). Such black holes are one of the interesting candidates 



for the sources of antiproton fiuxes recently observed in the BESS detector |14 



^For the case of n = 4, we do not find a consistent parameter region with the BESS experiment. 
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The dark matter PBHs play a role of the cold dark matter on the large scale structure 
formation. The scales of the fluctuations for PBH formation themselves are much smaller 
than the galactic scale and thus we cannot see any signals for such fluctuations in 5T /T 
measurements. However, the PBHs may be a source of gravitational waves. If the PBHs 
dominate dark matter of the present universe, some of them likely form binaries. Such 
binary black holes coalesce and produce significant gravitational waves [^] which may be 
observable in future detectors. 
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FIGURES 




0.2 0.4 0.6 0.8 

FIG. 1. The amplitude of density fluctuations (6^) and PBH masses. Here we take 
n = 4:, g = l,v = 10^''. The region on left hand side is irrelevant since a break scale exceeds 
our current horizon. The thick solid lines correspond to M-qh/Mq = 10^, 1, and 10~^, from left to 
right. The short dash - long dash lines show r^BH — 1 (top) and 0.1 (bottom). 
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0.2 0.4 0.6 0.8 

FIG. 2. The amplitude of density fluctuations {6b) and PBH masses. Here we take 
n = 3, g = 10~^,f = 10^^. The region on down right corner is irrelevant since a break scale 
is too small and does not cross horizon during inflation. The thick solid lines correspond to the 
break masses M^/Mq = 10'^^, 10'^'^, 3 x IQ-^*^, and IQ-^o^ from left to right. The short dash - 
long dash line shows Q{Mcva.p)h'^ ~ 2 x 10~^. The two thin solid lines show the region Q{Mi))h? ~ 1 
(top) and ~ 2 X 10~^ (bottom), respectively. 



14 



